Wealth distribution: To be or not to be a Gamma? 
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We review some aspects, especially those we can tackle analytically, of a minimal model of closed 
economy analogous to the kinetic theory model of ideal gases where the agents exchange wealth 
amongst themselves such that the total wealth is conserved, and each individual agent saves a 
fraction (0 < A < 1) of wealth before transaction. We are interested in the special case where the 
fraction A is constant for all the agents (global saving propensity) in the closed system. We show 
by moment calculations that the resulting wealth distribution cannot be the Gamma distribution 
that was conjectured in Phys. Rev. E 70, 016104 (2004). We also derive a form for the distribution 
at low wealth, which is a new result. 
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I. INTRODUCTION 

The distribution of wealth or income in society has 
been of great interest for many years. As first noticed by 
Pareto in the 1890's pQ, the wealth distribution seems to 
follow a "natural law" where the tail of the distribution 
is described by a power-law f(x) ~ x~( 1+a \ Away from 
the tail, the distribution is better described by a Gamma 
or Log-normal distribution known as Gibrat's law [2J. 
Considerable investigation with real data during the last 
ten years revealed that the power-law tail exhibits a re- 
markable spatial and temporal stability and the Pareto 
index a is found to have a value between 1 and 2 H] . 
Even after 110 years the origin of the power-law tail re- 
mained unexplained but recent interest of physicists and 
mathematicians in econophysics has led to a new insight 
into this problem (see Refs. [SHI]). 

Our general aim is to study a many-agent statistical 
model of closed economy (analogous to the kinetic the- 
ory model of ideal gases) [SHE], where N agents ex- 
change a quantity x, that may be defined as wealth. 
The states of agents are characterized by the wealth 
{xi}, i — 1, 2, . . . , N, and the total wealth W = J2i x i 
is conserved. The evolution of the system is then car- 
ried out according to a prescription, which defines the 
trading rule between agents. These many-agent statisti- 
cal models have N basic units {1,2,..., N}, interacting 
with each other through a pair-wise interaction charac- 
terized by a saving parameter A, with < A < 1. Wo 
define the equilibrium distribution of wealth f(x) as fol- 
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lows : f(x)dx is the probability that in the steady state 
of the system, a randomly chosen agent will be found to 
have wealth between x and x + dx. In these models, if 
A is equal for all the units, f(x) is fitted quite well by a 
Gamma-distribution fl4HT6l 
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This equilibrium distribution |T]) had been suggested by 
an analogy with the kinetic theory of gases in -D(A) di- 
mensions [T4Hl6j . 

In this paper we show by the method of moment 
calculations that the resulting wealth distribution can- 
not be the Gamma distribution that was conjectured in 
Ref . [T5] [T5] . We also derive the functional form of an 
upper bound on f(x) at very small x. 



II. MANY-AGENT MODEL OF A CLOSED 
ECONOMY 

We study many-agent statistical models of closed econ- 
omy (analogous to the kinetic theory model of ideal 
gases), where N agents exchange wealth x. The states 
of agents are characterized by the wealth {x{\, i = 
1, 2, . . . , N, and the total wealth W = Xi is conserved. 
The evolution of the system is then carried out accord- 
ing to a prescription, which defines the trading rule be- 
tween agents. At every time step two agents i and j 
are extracted randomly and an amount of wealth Ax is 
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exchanged between them, 



c' = Xj — Aa 



Ax. 



(3) 



It can be noticed that in this way, the quantity x is con- 
served during the single transactions: x' { + x'j = Xi + Xj 
(see Fig. [TJ) , where x\ and x'j are the agent wealths after 
the transaction has taken place. Several simple models 
dealing with different transaction rules have been studied 
(see the reviews [5j El [17] and references therein) . Here 
we will present a few examples. 
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A. Basic model without saving: Boltzmann 
distribution 



In the first version of the model, the wealth difference 
Air is assumed to have a constant value [5], 

Ax = Ax . (4) 

This rule, together with the constraint that transactions 
can take place only if x\ > and x'j > 0, provides a 
Boltzmann distribution, see the curve for A = in Fig. [2] 
Alternatively, Ax can be a random fraction of the wealth 
of one of the two agents, 

Ax = exi or Aa; = exj , (5) 

where e is a random number uniformly distributed be- 
tween and 1. A trading rule based on the random re- 
distribution of the sum of the wealths of the two agents 
had been introduced by Dragulescu and Yakovenko [TT] . 

x^ ^{xi ~h Xj) , 

x'j = (1 - e)(xi + xj) . (6) 

Equations Q are easily shown to correspond to the trad- 
ing rule ([3]), with 

Aa; = (1 — e)xi — exj . (7) 

All the versions of the basic model lead to an equilibrium 
Boltzmann distribution, given by 

/(x) = _l exp (_JLj , (8) 

where the effective temperature of the system is just the 
average wealth (x) [5] [TT]. The result mty is found to 
be robust; it is largely independent of various factors. 
Namely, it is obtained for the various forms of Aa; men- 
tioned above, for a pair-wise as well as multi-agent inter- 
actions, for arbitrary initial conditions |12) . and finally, 
for random or consecutive extraction of the interacting 
agents. For the trading rule Q one can show the conver- 
gence towards the Boltzmann distribution through dif- 
ferent methods: Boltzmann equation, entropy maximiza- 
tion, distributional equation, etc. 



FIG. 1. Analogy of the minimal economic model with a clas- 
sical isolated system of ideal gas, where the particles are ran- 
domly undergoing "Elastic" collisions, and exchanging kinetic 
energy. In the closed economy, the economic agents randomly 
trade with each other according to some rule and exchange 
wealth. 



B. Model with global saving propensity A 

A step toward generalizing the basic model and making 
it more realistic, is the introduction of a saving criterion 
regulating the trading dynamics. This can be practically 
achieved by defining a saving propensity < A < 1, 
which represents the fraction of wealth which is saved - 
and not reshuffled - during a transaction. The dynamics 
of the model is as follows [HI [TH] : 

x[ = Xxi + e(l — X)(xi + Xj) , 

x'j = Xxj + (1 - e)(l - X)(xi + Xj) , (9) 

corresponding to a Aa; in Eq. ^ given by 

Ax=(l-X)[(l-e)xi-exj]. (10) 

This model leads to a qualitatively different equilibrium 
distribution. In particular, it has a mode x rn > and 
a zero limit for small x, see Fig. [2] Later we will derive 
a form for an upper bound on fix) at low range. The 
functional form of such a distribution was conjectured to 
be a T-distribution, as given by Eq. (JlJ on the basis of 
an analogy with the kinetic theory of gases. Indeed, it is 
easy to show, starting from the Maxwell-Boltzmann dis- 
tribution for the particle velocity in a I? dimensional gas, 
that the equilibrium kinetic energy distribution coincides 
with the Gamma-distribution ([!]) with n = y . This con- 
jecture is remarkably consistent with the fitting provided 
to numerical data [MHTB] . In the following section we will 
show by two different approaches that the conjecture ([!]) 
cannot be the actual equilibrium distribution. 
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FIG. 2. Probability density for wealth x. The curve for A = 
is the Boltzmann function f(x) — (x)^ 1 exp(— x/{x)) for the 
basic model of Sec. |II A| The other curves correspond to a 
global saving propensity A > (see Sec. |II B[ ). 



III. ANALYTICAL RESULTS FOR MODEL 
WITH SAVING PROPENSITY A 

A. Fixed-point distribution 

Let X be a random variable which stands for the 
wealth of one agent, at equilibrium and in the limit of 
an infinite number of agents, we can say from Eq. ^ 
that the law of X f, is a fixed-point distribution of the 
equation 



X = XX 1 +e(l-X)(X 1 +X 2 ), 



(11) 



where = means identity in distribution and one assumes 
that the random variables X\ , X 2 and X have the same 
probability law, while the variables Xi,X 2 and e are 
stochastically independent. It seems difficult to find the 
distribution of X, however, one can compute the mo- 



ments of /. Indeed with ( |11[ ), one can write immediately 
Vm e N, (X m ) = ((XX! + e(l - X)(X X + X 2 )) m ) , (12) 



and by developing ( 12 ) one can find the recursive relation 

v x m ~ k (l - X) k A (k 
k J k + 1 ^ \p 



(X r ' 



k=0 



p=0 



(x m - p ) (X p ) 



(13) 

Using (13 1 with initial conditions (X°) = 1 (normal- 
ization) and (X 1 ^ = 1 (without loss of generality), we 
obtain 



<* 2 > 

(X*) 
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1 + 2A ' 
3(A + 2) 
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(14) 
(15) 
(16) 



Now let us compare theses moments with conjecture Q's 
moments. Setting (x) = 1 in Eq. ([!]) it is easy to show 

^ = (n + fc-l)(n + fc-2)...(n+l) ^ 



Writing (17) for k — 2,3,4 and choosing n as in ^ we 
find 



<* 3 > = 



1 



A + 2 



3(A + 2) 



(18) 
(19) 



n 1 + 2A 

(n + 2)(n+ 1) _ 

n 2 ~ (1 + 2A) 2 ' 

(n + 3)(n + 2)(n+l) _ 3(A + 2)(4-A) 

n 3 ~ (1 + 2A) 3 ■ [ ' 



The fourths moments (eqs.(16l and (p0|) are different so 



the conjecture that the Gamma distribution is an equi- 
librium solution of this model is wrong. Nevertheless the 
first three moments coincide exactly which shows that 
the Gamma-distribution is strangely a very good approx- 
imation. Moreover the deviation in the fourth moment is 
very small (see Fig. |3j which shows that the two curves 
can hardly be distinguished by the naked eye). Finding 
a function that would coincide to higher moments is still 
an open challenge. These results are consistent with the 
ones found by Repetowicz et al. [T5] which will be pre- 
sented in the following section. 



B. Laplace transform analysis 

In this section we will confirm the previous result with 
a different approach based on the Boltzmann equation 
and along the lines of Bassetti et al. [20 . Given a fixed 
number of N agents in a system, which are allowed to 
trade, the interaction rules describe a stochastic process 
of the vector variable (xi(t), . . . , xn(t)) in discrete time 
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FIG. 3. Exact fourth moment eq.(16l and Gamma distribu- 



tion fourth moment eq.(20l against A. The inset shows the 



relative difference between exact fourth moment eq.(I6l and 



Gamma distribution fourth moment eq.(20l against A 



r. Processes of this type are thoroughly studied e.g. in 
the context of kinetic theory of ideal gases. Indeed, if the 
variables Xi are interpreted as energies corresponding to 
the i-th particle, one can map the process to the mean- 
field limit of the Maxwell model of particles undergo- 
ing random elastic collisions. The full information about 
the process in time r is contained in the ./V-particle joint 
probability distribution Pn{t, Xi,X2, ■ ■ ■ ,%n)- However, 
one can write a kinetic equation for one-marginal distri- 
bution function 

Pi(t,x) = J Pn(t,x,x 2 , ■ ■ ■ , XN)dx 2 ■ ■ ■ dxN, 

involving only one- and two-particle distribution func- 
tions 

Pi(t+1,x)-Pi(t,x) = J P 2 {t,x,, Xj ) 

(s(x — Xx — (1 — X)c(xi + X j)) 
+ S(x — Ax — (1 — A)(l — e)(xi + Xj))jdxidxj 

-2P 1 (r,x) 

which may be continued to give eventually an infinite 
hierarchy of equations of BBGKY (Born, Bogoliubov, 
Green, Kirkwood, Yvon) type [3T]. The standard ap- 
proximation, which neglects the correlations between the 
wealth of the agents induced by the trade gives the fac- 
torization 

P 2 (T, Xi ,Xj) = Pi(T,Xi)Pi(T,Xj), 

which implies a closure of the hierarchy at the low- 
est level. In fact, this approximation becomes exact 



in the thermodynamic limit (N oo). Therefore, the 
one-particle distribution function bears all information. 
Rescaling the time as t = ^ in the thermodynamic limit 
N — > oo, one obtains for the one-particle distribution 
function f(t,x) the Boltzmann-type kinetic equation 



df(t,x) 

at 



f(t,Xi)f(t,Xj) 

(j)(x — Xx — (1 — X)e(xi + X j)) 

+ S(x — Xx — (1 — A)(l — e){xi + Xj))j dxi dx^j 

-f(t,x). (21) 

This equation can be written (see Matthes et al. [20]) as 

df(t,x) 



dt 



Q(f,f), 



where Q is a collision operator. A collision operator is 
bilinear and satisfies, for all smooth functions (j)( x ) 



Q(fJ)<t>(x)dx 



OO rOO 



(4>{x' i ) + <l>{x' i )-4 > {x i )-cl>{x j )) 



f(x i )f(x j )dx l dx j ) 1 (22) 



where x\ and x'j are the post-trade wealth. With this 
property the equation can be written in the weak form, 
for all smooth functions d>(x) 



dt 



f(t, x)<p{x)dx 



OO f>00 



f{xi)f(xj)d Xl d Xj ). (23) 



e sx gives 



It is very useful because the choice <j){x) 
(after some calculations) the Boltzmann equation for the 
Laplace transform / of / 



f(t,s) 



f(t,(X+(l-X)e)s)f(t,(l-X)es) 



df(t,a) 
dt 



+ f(t, (1 - A)(l - e)s)f(t, 1 - (1 - X)es)). (24) 

For the steady state, and if e is drawn randomly from a 
uniform distribution, the previous equation reduces to 



a/00 



1 



1 - A 



(1-A) S 



/(As + y)f{y)dy, (25) 



which coincides with results of |19) . The Taylor expan- 
sion of /(s) can be derived by substituting the expan- 



sion f(s) = EZv(- 1 Y m pS P in (25). Since f(-s) is 



the moment-generating function we have (x k ) = mu ■ k\. 
With this method Repetowicz et al. [H] obtained the 
recursive formula 



with 



nip = m„in 



Cf\\) 
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1 - A 



(p) _ {i-\f-^-{ q + i)c { q p) 



°<7 + l 



p-q 



(26) 



x/(l-X) 




x/X x 



~ {p) = [i-xy 
p + i 



(27) 



Now with this formula one can obtain the first four mo- 
ments and they match the ones found in the previous 
section eqs. (14-16), which confirms that the Gamma- 
distribution is not the stationary distribution. 



FIG. 4. Region of integration 



However, one can simplify the equation, by doing the in- 
tegral over e. Then the <5-function will contribute only if 
we have the following constraints 



< Xi < x/X, 



(31) 



IV. UPPER BOUND FORM AT LOW WEALTH 
RANGE 



1-A 



(32) 



From equation (111 



X = XX t + e(l — X)(Xi + Xj), 

we have for all x > 

F[X <x]= FiXX, + e(l - X)(X l + Xj) < x], (28) 

where P[.] means the probability of the event inside the 
brackets. If the number of agents in the market is 
large, the distributions of different agents are indepen- 
dent. Then 

/>oo poo 

dxf(x) = / dxif{xi) / dxjfixj) 
Jo Jo 

x / deQ[x — Xxi + e(l — X)(xi + xj)] , 
Jo 

(29) 

where <d is the Heaviside step function. Taking the 
derivative with respect to x in both sides, we have 



oo 



3G 



f(x) = j dx l f{x i ) J dxjf(xj) 

[ deS[ 
Jo 



JO 

I <i<<)\.r- Xxi-e(l- X)(xi + Xj)]. (30) 



This equation is an integral equation for f(x). As men- 
tioned earlier, we are not able to solve it in closed form. 



< x-j 



(33) 



The range defined by these constraints is shown in figure 
[4] In this range, the derivative of the argument of the 
delta function with respect to e is just (x t + Xj)(l — A). 
And, hence we get 



/(*) = 



1 



This immediately gives 



/A /-oo ^ 

dxif{xi) / _ dxjf(xj) — - — 

Jmax(^J,0) Xi+Xj 



x/X 



where 



fix) < C / f{xi)dxi 
Jo 



C= / dxjffa)—. 



(34) 
(35) 

(36) 



We assume that / decays fast enough near so that 



the integral in (36) is well defined. Now (35) may be 



rewritten by rescaling the variable, as 



f{Xx) < C / dxif{xi). 
Jo 



(37) 



We now use the observation that for A > the numer- 
ically determined f(x) is a continuous function with a 
single maximum, say at xq (see Fig. ^J. Then for all 
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x < xq, the integrand (37 1 takes its maximum value at 



the right extreme point, i.e. when Xi = x. This then 
gives us 

f(Xx) < Cxf(x), for x < x . (38) 
Iterating this equation, we get 

f(\ r x) < C r \ r{r - 1)/2 x r f{x). (39) 
We can set x = x in the above equation, giving 

f(X r x ) < C r X r ^/ 2 x r f(x ). (40) 



Then taking r 
get 



log x and rescaling the variables, we 



f(x)^0(x a exp[-f3(\ogx) 2 }) 



(41) 



as x — > 0, where a and /?(> 0) are two constants depen- 
dent on A. The Gamma-distribution decays slower than 
the rhs in (41) when x — > 0. The expression (41) gives 



an upper bound form at low wealth range and confirms 
again that the distribution of the global saving propen- 
sity model is not a Gamma-distribution. 



V. DISCUSSION AND OUTLOOK 

We have used different approaches to show that the 
correct form of the wealth distribution cannot be the 
Gamma distribution. We have also derived an analyt- 
ical form of an upper bound at low wealth range see Eq. 
(41 ). This is an analytically calculated upper bound but 



the closed form of the solution to Eq. ^ is still an open 
question. 

As a further generalization, the agents could be as- 
signed different saving propensities Xi [131 Q3B |2"2"H2"5] . In 
particular, uniformly distributed Xi in the interval (0, 1) 
have been studied numerically in Refs. [HI [22] • This 
model is described by the trading rule 

x\ = XiX. t + e[(l - Xt)xi + (1 - Xj)Xj] , 

x] = XjXj + (1 - e)[(l - Xi)xi + (1 - X 3 ) Xj ] , (42) 

or, equivalently, by a Aa; (as defined in Eq. (|3|) given by 

Ax = (1 - e)(l - Xi)xi - e(l - X 3 )x 3 . (43) 

One of the main features of this model, which is sup- 
ported by theoretical considerations [HI [23l [26] , is that 
the wealth distribution exhibits a robust power-law at 
large values of x, 



(44) 



with a Pareto exponent a = 1 largely independent of the 
details of the A-distribution. 

ACKNOWLEDGMENTS 

The authors are grateful to D. Dhar for critical read- 
ing of the manuscript and inputs in finding the functional 
form of the wealth distribution at low range. The authors 
also acknowledge F. Abergel, N. Millot, M. Patriarca, M. 
Politi and A. Chatterjee for critical discussions or com- 
ments, and B.K. Chakrabarti for pointing out reference 
[7J. AC is grateful to Department of Theoretical Physics, 
TIFR for the kind hospitality where part of the work was 
initiated. 



[1] V. Pareto, Cours d'economie politique, F. Rouge, Lau- 
sanne (1897). 

[2] R. Gibrat, Les inegalites economiques, Sirey, Paris 
(1931). 

[3] A. Dragulescu, V.M. Yakovenko, Exponential and power- 
law probability distributions of wealth and income in the 
United Kingdom and the United States, Physica A 299 
(2001) 213. 

[4] H. Aoyama, W. Souma, Y. Fujiwara, Growth and fluctu- 
ations of personal and company's income, Physica A 324 
(2003) 352. 

[5] V. Yakovenko, J. Rosser, Colloquium: Statistical me- 
chanics of money, wealth, and income, Rev. Mod. Phys. 
81 (2009) 1703. 

[6] A. Chatterjee, B.K. Chakrabarti, Kinetic Exchange Mod- 
els for Income and Wealth Distributions, Eur. Phys. J. B 
60 (2007) 135. 

[7] A.S. Chakrabarti and B.K. Chakrabarti, Statistical the- 
ories of income and wealth distribution, Economics E- 
Journal 4 (2010) 2010-4. 

[8] E. Bennati, La simulazione statistica nell'analisi della 
distribuzione del reddito: modelli realistici e metodo di 



Monte Carlo, ETS Editrice, Pisa (1988); E. Bennati, 
Un metodo di simulazione statistica nell'analisi della dis- 
tribuzione del reddito, Rivista Internazionale di Scienze 
Economiche e Commerciali 35 (1988) 735; E. Bennati, II 
metodo Monte Carlo nell'analisi economica, Rassegna di 
lavori dell'ISCO X (1993) 31. 
[9] J. Angle, The Surplus Theory of Social Stratification and 
the Size Distribution of Personal Wealth, Social Forces 65 
(1986) 293; J. Angle, The inequality process as a wealth 
maximizing process, Physica A 367 (2006) 388. 
[10] S. Ispolatov, P.L. Krapivsky, S. Redner, Wealth distribu- 
tions in asset exchange models, Eur. Phys. J. B 2 (1998) 
267. 

[11] A. Dragulescu, V.M. Yakovenko, Statistical mechanics of 

money, Eur. Phys. J. B 17 (2000) 723. 
[12] A. Chakraborti, B.K. Chakrabarti, Statistical mechanics 

of money: How saving propensity affects its distribution, 

Eur. Phys. J. B 17 (2000) 167. 
[13] A. Chatterjee, B.K. Chakrabarti, S. S. Manna, Money in 

gas-like markets: Gibbs and Pareto laws, Physica Scripta 

T 106 (2003) 367. 
[14] A. Chakraborti, M. Patriarca, Gamma-distribution and 



7 



Wealth inequality, Pramana 71 (2008) 233. 
[15] M. Patriarca, A. Chakraborti, K. Kaski, Statistical 

model with a standard gamma distribution, Phys. Rev. 

E 70 (2004) 016104. 
[16] M. Patriarca, A. Chakraborti, K. Kaski, Gibbs versus 

non-Gibbs distributions in money dynamics, Physica A 

340 (2004) 334. 

[17] M. Patriarca, E. Heinsalu, A. Chakraborti, Basic kinetic 

wealth-exchange models: common features and open 

problems, Eur. Phys. J. B 73 (2010) 145. 
[18] A. Chakraborti, Distributions of money in model markets 

of economy, Int. J. Mod. Phys. C 13 (2002) 1315. 
[19] P. Repetowicz, S. Hutzler, P. Richmond, Dynamics of 

money and income distributions, Physica A 356 (2005) 

641. 

[20] D. Matthes, G. Toscani, On steady distributions of 
kinetic models of conservative economies, J. Stat. 
Phys. 130 (2008) 1087; F. Bassetti, G. Toscani, 
Explicit equilibria in a kinetic model of gambling, 
http://arxiv.org/abs/1002.3689l 



[21] M. Plischke, B. Bergersen, Equilibirum statistical me- 
chanics, World scientific (1989) 137. 

[22] A. Chatterjee, B.K. Chakrabarti, S. S. Manna, Pareto 
law in a kinetic model of market with random saving 
propensity, Physica A 335 (2004) 155. 

[23] A. Chatterjee, B.K. Chakrabarti, R. Stinchcombe, Mas- 
ter equation for a kinetic model of trading market and 
its analytic solution, Phys. Rev. E 72 (2005), 026126. 

[24] M. Patriarca, A. Chakraborti, K. Kaski, G. Germano, 
Kinetic theory models for the distribution of wealth: 
Power law from overlap of exponentials, in: A. Chat- 
terjee, S.Yarlagadda, B. K. Chakrabarti (Eds.), Econo- 
physics of Wealth Distributions, Springer, (2005) 93. 

[25] A. Chakraborti, M. Patriarca, Variational principle for 
the Pareto power law, Phys. Rev. Lett. 103 (2009) 
228701. 

[26] P.K. Mohanty, Generic features of the wealth distribution 
in an ideal-gas-like market, Phys. Rev. E 74(1) (2006) 
011117; U. Basu, P.K. Mohanty, Modeling wealth in 
growing markets, Eur. Phys. J. B 65 (2008) 585. 



